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Abstract 

We investigate electronically excited atoms in a magnetic guide. It turns out that the Hamilto- 
nian describing this system possesses a wealth of both unitary as well as antiunitary symmetries 
that constitute an uncommon extensive symmetry group. One consequence is the two- fold de- 
generacy of any energy level. The spectral properties are investigated for a wide range of field 
gradients and the spatial distributions of the spin polarization are analyzed. Wave lengths, oscil- 
lator strengths and selection rules are provided for the corresponding electromagnetic transitions. 
The effects due to an additional homogeneous bias field constituting a Ioffe-Pritchard trap are 
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I. INTRODUCTION 



External fields are nowadays widely used to control the motion of atoms including their 
cooling and trapping as well as the preparation of their internal states. Optical lattices and 
atom chips are two major examples of devices that allow to deal with atomic ensembles 
but also possess the perspective of manipulating single atoms for the purpose of quantum 
information processing. To this end it is indispensable to understand the structure and 
behaviour of (excited) individual atoms in traps. In the case of the atom chip (see ref . |l( and 
references therein) tight magnetic traps on the micrometer scale can be created exhibiting 
large field gradients which are not accessible in the case of macrosopic traps. Highly excited 
Rydberg atoms therefore start to 'feel' the variation i.e. the inhomogeneity of the magnetic 
field across the extension of their wave functions. This naturally leads to the question: How 
do inhomogeneous magnetic field configurations alter the electronic structure of excited 
atoms ? 

During the past decades many thorough investigations have been performed on the be- 
haviour and properties of excited (Rydberg-) atoms in homogeneous magnetic fields (see the 
books and reviews DDDDD . Indeed investigations on atoms in strong magnetic fields 
provided major contributions to a variety of different research areas such as semiclassics 
of nonintegrable systems, 'quantum chaos', nonlinear dynamics, astrophysics of magnetized 
stars and it elucidated and significantly advanced our understanding of magnetized struc- 
tures in general. 

In contrast to the case of a homogeneous magnetic field there exist no studies on the 
electronic structure of atoms in the presence of inhomogeneous external fields: all investiga- 
tions in the literature on the behaviour of ultracold atoms in inhomogeneous fields typically 
treat the atom as a point particle whose magnetic moment couples either adiabatically or 
nonadiabatically [?], 0, 10 L 11, ll^ to the external field. This holds with the exception of 



two very recent works |l3. Il7| that consider the electronic structure of atoms with a single 
active electron subject to a three-dimensional quadrupole field. A variety of interesting new 
phenomena have been observed there. The symmetries of this system cause each energy 
level to be degenerate in the presence of the field. Furthermore the intimate coupling of 
the spin and spatial degrees of freedom leads to a complex spatial distribution of the spin 
polarization of individual electronic states. A remarkable property of the electronic states 
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in the 3D-quadrupole trap is the fact that they possess a magnetic field-induced permanent 
electric dipole moment whose size strongly varies with the Rydberg state considered. Besides 
the 3D-quadrupole field there is another generic inhomogeneous magnetic field configuration 
which is employed to trap atoms in particular on the atom chip Q . This is the so-called side 
guide which is created by superimposing the magnetic field of a current carrying wire with 
a homogeneous bias field oriented perpendicular to the wire. The resulting magnetic guide 
can be augmented to a Ioffe-Pritchard type 2D-trap by applying an additional homogeneous 
bias field parallel to the wire. It is exactly this configuration which is studied in the present 
work i.e. we investigate the structure and properties of electronically excited atoms in a 
magnetic guide. According to the effects obtained for atoms in a 3D-quadrupole trap in 
Refs. [l^, [itJ we expect also the atoms in a side guide to exhibit interesting new features. 

The paper is organized as follows: In Sec. |H] we introduce the field configuration gen- 
erated by a so-called side guide. We specify our approach which is particularly suited for 
ultra-cold atoms with a single active electron and derive the corresponding Hamiltonian. 
This Hamiltonian exhibits a wealth of both unitary and anti-unitary symmetries and consti- 
tutes an uncommon large symmetry group which is analyzed in Sec. IIHI In particular these 
symmetries lead to a two-fold degeneracy of any energy level, similar to the case of an atom 
in a 3D-quadrupole trap. A discussion of an arbitrary spin-|-systems in a field configuration 
obeying certain symmetries are discussed. Section |IV| contains a discussion of the properties 
of the symmetry-adapted electronic states. In Sec. El the latter are studied in case an addi- 
tional homogeneous (Ioffe-)field is applied. The numerical methods being employed in order 
to solve the stationary Schrodinger equation are briefly outlined in Sec. IVII A discussion of 
our results is provided in Sec. IVHI We analyze the spectra for a wide range of gradients. 
Furthermore we explore properties of the electronic spin such as spin expectation values and 
distributions of the spin polarization. Selection rules and dipole strengths of electric dipole 
transitions are calculated. We close with a discussion of the electronic structure in case a 
homogeneous magnetic field is applied in addition to the field of the magnetic guide. Sec. 
IVII II contains the summary and outlook. 
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II. THE FIELD CONFIGURATION AND THE HAMILTON! AN 



Alkali atoms are used throughout many experiments in ultra cold atomic physics. Besides 
a single active electron they possess a closed shell core and the total electronic spin is 
therefore exclusively carried by the outer electron. We assume the motion of this valence 
electron to take place in the Coulomb potential of a single positive point charge. Since the 
focus of this work is to understand fundamental features of electronically excited atoms in a 
certain inhomogeneous magnetic field we do not account for quantum defects which would 
require the consideration of core-electron scattering processes. We also neglect relativistic 
effects such as spin-orbit and hyperfine coupling. Both interaction possess a r _3 -dependence 
with r being the distance between the outer electron and the nucleus. For (highly) excited 
states their contributions can safely be neglected or, if necessary accounted for by means 
of perturbation theory. Since we focus on ultra cold atoms effects of the center of mass 
(cm.) motion on the electronic motion are neglected here. Specifically we assume an 
infinitely heavy core (cm.) located at the minimum of the magnetic field. Employing the 
above approximations the Hamiltonian describing the motion of the valence electron in the 
presence of an external magnetic field reads 
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The magnetic field is introduced via the minimal coupling including the vector potential 
thereby providing the kinetic energy in the presence of the field. The third term represents 
the coupling between the spin of the electron and the external field. A common configuration 

n 

for the manipulation of neutral atoms is the so-called magnetic side guide jlj . This particular 
setup is generated by a current carrying wire whose 'circular' magnetic field is superimposed 
by an external homogeneous bias field perpendicular to the current flow. As a result the field 
vanishes along a line parallel to the wire at a distance po — being completely determined 
by the current I and the homogeneous magnetic field strength B. The Taylor expansion of 
the field around p yields 
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These are the quadrupolar, hexapolar and octopolar expansion terms of the field. Here 
we restrict ourselves to the linear term which should provide a good approximation of the 
magnetic field configuration as long as p 3> 1. Thus we obtain the expression 
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(3) 



Here b is the magnetic field gradient determining the linear growth of the field with increasing 
distance from the line of zero field. Figure Q shows two vectorial plots along cuts through 
the field. The cut through the x — y-plane reveals the quadrupolar shape of the field of the 
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FIG. 1: Vectorial plots of the magnetic field ©• a: Intersection for z = 0. The quadrupolar shape 
of the field is clearly recognized, b: Intersection for x = revealing the translational invariance 
with respect to the z-coordinate. 



side guide whose translational invariance along the ,2-axis can be easily observed in figure 
\I]p. A corresponding vector potential in the Coulomb gauge is given by 



A = b 







(4) 



nserting the expressions (jSJ) and (0J) into the Hamiltonian |Q thereby adopting atomic units 
yields 



1 1 
H = -A 7 = 

2 y/ x 2 + y 2 + z 
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(5) 



The first two terms of (J^J) represent the non-relativistic hydrogen atom. The third term 
which is linear with respect to b replaces the angular Zeeman term 2(| which would occur 
in a homogeneous field. Here the spatial coordinates in x and y couple with the momentum 
in z direction. The successive diamagnetic term oc b 2 represents an oscillator coupling term 
confining the electronic motion in the x and y direction except for the axis exit channels. 
This is reminiscent but also very different to the situation in a homogeneous field, where 
the diamagnetic interactions in the x and y direction separate and represent pure harmonic 
oscillators. Finally the fifth term represents the coupling of the electronic spin to the spatial 
coordinates arises from the interaction of its magnetic moment with the field. We here 
encounter a linear dependence on the spatial coordinates and the gradient b. This term 
prevents the factorization of the motions in coordinate space and spin space. Finally one 
should note that the only explicit dependence on the coordinate z is due to the Coulomb 
term. Without this rotationally invariant interaction the system would be invariant under 
translations with respect to the ^-coordinate. 

Performing the canonical scaling transformation x = b^x and p = b~sp the Hamiltonian 
(jSJ becomes 



1 -2 % 1 2 2 1 



H = b *H=-p + xyp z + - x y + - (xa x - ya y ) (6) 

1 v x + y + z 

with Z = and where we have for simplicity omitted the bar on top of the phase space 
variables. This shows us that employing a scaled energy (scaled Hamiltonian) the only free 
parameter is the scaled Coulomb coupling strength Z that depends on the field gradient. 
The scaled Hamiltonian describes the motion of an electron in the Coulomb-field of a charge 
Z and the field with gradient 1. If b — > oo the Coulomb term vanishes since Z —>■ 0. In this 
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limit the energy level spacing is expected to scale according to &3 . 



III. SYMMETRIES AND DEGENERACIES IN SPIN-±-SYSTEMS 

In this section we analyze the structure of the Hamiltonian (j3J) in detail. After studying 
its symmetries we discuss how these symmetries affect the excitation spectrum. As a result 
of a tedious and elaborate analysis of the Hamiltonian (jSJ) we found 15 distinct symmetry 
operations leaving it invariant. A complete list is provided in table |U Each symmetry is 
composed of a number of elementary operations which are shown in table [H] All symmetry 
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TABLE I: Symmetry operations of the Hamiltonian ©. Top part: unitary symmetries. Bottom 
part: anti-unitary symmetries. 
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TABLE II: Set of discrete operations out of which all symmetry operations of the Hamiltonian © 
can be composed. Note that S\ and S2 are given in a basis where a z is diagonal. 

operations are either unitary or anti-unitary. The anti-unitary ones involve the conventional 
time reversal operator T. In spite of its simplicity our system therefore possesses a wealth 
of symmetry properties. The algebra of the underlying symmetry group possesses a com- 
plicated structure some features of which will be discussed in the following. The operators 
Sjj., £ y and T> z generate a sub-group obeying the algebra [£j,£j] = 2iejjfc£& reminiscent of 
angular momentum operators. We have E? = 1. Interestingly these quantities act on both 
real and spin space. A deeper look into the representation theory of our group reveals a 
two-fold degeneracy of any energy level similar to those we encountered in our investigations 
of atoms in a three-dimensional quadrupole trap 
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Alternatively this degeneracy can also be established as follows. The operations S z and 
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Ta x P z obey {£ z , Tcr^-P^} = 0. Let \E,tt) be an energy eigenstate and at the same time an 
eigenstate of S z with 

Z z \E,ir) = tt\E,?t) (7) 
and 7r = ±1. Employing the above anti-commutator one obtains 

E z Ta x P z \E,7v) = -Ta x P z E z \E, vr) = -itT<j x P z \E, tt) (8) 

n 

The state Ta x P z \E, tt) can be identified with \E, — tt). Hence, as long as tt ^ [27] there is 
always an orthogonal pair of states possessing the same energy namely \E, tt) and \E, —tt). 
We have to emphasize that there occur no further degeneracies in the system. In principle 
one could think of performing the above calculation repeatedly but now substituting Ta x P z 
with any operator listed in table |l] which anti-commutes with S z . It turns out that all of the 
resulting states generated by this scheme are either superpositions of \E, tt) and \E, —tt) or 
differ only by a phase factor from one of these states. 

Out of the 15 symmetry operations one can choose several sets of commuting operators. 
For the following investigation we choose the set H, £ z , P y P z I xy S2- The combination of S z 
and PyP z I xy S 2 leads to the additional commuting operator P x P z I xy S%. We have found the 
properties: 

{PyPzIxyS2) = (PxPz^xyS^) = (9) 

(S z ) 2 = {P y PJ xy S 2 f = {P x PJ xy Slf = 1. (10) 

For completeness we provide here the general embedding of the above-derived degeneracies 
due to symmetries. Let us assume we have a general spin-|-systems with the following 
accompanying properties: 

1. There are two operators A and B commuting with the underlying Hamiltonian: 
[H,A] = [H,B] = 0. 

2. A and B anti-commute: {A, B} = 0. 

3. A is a Hermitian operator. B is an ( ant i-) unitary operator which can be written as a 
product B = RS where R and S exclusively act on the real space and the spin space, 
respectively. 
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4. The operator S is trace-less: Tr S = 0. 

If these conditions are fullfilled any state is doubly degenerate. This is seen as follows. 
Property 4 immediately leads to Tr B = 0. Hence, we find the nonzero eigenvalues of B to 
appear pairwise with opposite signs. If now \E, b) is an eigenstate of B and at the same 
time an energy-eigenstate property 2 implies that 

BA \E,b) = -AB \E, b) = -bA \E,b) = -b \E, -b) . (11) 

Hence, \E, b) and A \ E, b) = \E, —b) are two degenerate energy-eigenstates of the system. 

In the present case the two anti-commuting operators are S z and Ta x P z . For the case 
of an atom in a three-dimensional quadrupole field we have A = J z and B = Ta x P z . In 
a homogeneous magnetic field the remaining symmetries constitute an Abelian symmetry 
group leading to exclusively one dimensional irreducible representations i.e. no degeneracies 
occur. Finally we remark that the reader can find in ref. ^| a discussion of degeneracies 
in spin-i-systems based on the properties of time-reversal operators. 

IV. S 2 -, P y PJ xy S 2 - AND Tct^-EIGENSTATES 

The operator P y P z I xy S 2 obeys the eigenvalue relation 

PyPJ xy S 2 \K) = k\k) . (12) 

Since 

\K) = {P y P z I xy S 2 f\K) = ^\K) (13) 

the eigenvalue k can adopt the four values ±1 and ±z. The reader should note that P y P z I xy S 2 
is a unitary but non-Hermitian operator. We therefore encounter complex eigenvalues. If 
we apply T, z to the states |ac) we find by exploiting equation © 

S z | K = ±i) = \ K = ±i) (14) 
S 2 \k = ±1) = - \k = ±1) . (15) 

By using the relation (Tcr x P z ) (P y P z I xy S 2 ) — i (P y P z I xy S 2 ) (Ta x P z ) = one finds the degen- 
erate pairs of states in the P^Pji^S^-subspaces: \E, +1),\E, —i) and \E, — 1),|P, Since 
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non-Hermitian operators do not represent physical observables only the quantum number n 
should be of direct relevance for the experimental observation. 

We now derive the expectation value of an observable Y in an eigenstate of E a . Assume 
we have {Y, T, z } = and hence 

(E,n\Yi: z \E,tt) = -{E,n\Z z Y\E,7i} (16) 
n(E,n\Y \E,n) = -tt {E,n\Y \E,tt) . (17) 

This immediately leads to the result 

(E,tt\Y\E,it) = 0. (18) 

The same arguments hold for an observable Z obeying {Z, P y P z I xy S2} = in which case we 
obtain 

(E,k\Z\E,k) = 0. (19) 

In the preceding section we showed the degeneracy of the states \E, n) and Ta x P z \ E,n). 
By superimposing these two states eigenstates of the operator Ta x P z can be constructed: 

\E, ±) Tu * P > = -j= [\E, tt) ± Ta x P z \E, tt)] . (20) 

The corresponding eigenvalue relation is 

Ta x P z \E,±) = ±\E,±). (21) 



V. ADDITIONAL HOMOGENEOUS FIELD IN z-DIRECTION (IOFFE FIELD) 

The application of an additional homogeneous magnetic field along the ^-direction (Ioffc 
field) has a dramatic impact on the properties of the system. In particular the symmetry 
properties are affected. The Hamiltonian becomes 

Hi = ~\^y,z - , 1 = + bxyp z + h - x 2 y 2 + \ (xa x - ya y ) 
* \j x 1 + y 2 + z 2 * * 

(xp v - yp x ) + ^-( x2 + y 2 ) + Y az (22) 

with B[ being the field strength of the Ioffe field. Since both the 2D-quadrupole (due to 
the side guide) and the magnetic field are perpendicular to each other the homogeneous 
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field terms can simply be added to the Hamiltonian (J3J). We find the well known Zeeman 
as well as the diamagnetic oscillator term. The coupling of the spin to the Ioffe field leads 
to a term being proportional to a z . The symmetries of Hj are listed in table ITTTl Due to 
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TP x a z 


TP z I xy S2 


PPx Py Pz Ixy &2 


TPy 



TABLE III: Symmetries of the Hamiltonian (|22jl. i.e. side guide with Ioffe field. Top line: unitary 
symmetries. Bottom line: anti-unitary symmetries. 

the presence of the additional homogeneous field numerous symmetries are lost (see table 
for comparison). The remaining operations form a non-Abelian algebra. In contrast to the 
group operations listed in table H] there are no two anti-commuting operators. Hence it is 
not possible to construct pairs of degenerate energy eigenstates as discussed above. Thus, 
applying the Ioffe field lifts the degeneracies occuring in the absence of it. Even with a 
finite Ioffe field the operations S 2 , P y P z I xy S% and P x P z l xy S\ together with Hi form a set of 
commuting operators. 



VI. NUMERICAL TREATMENT 



In order to obtain many eigenvalues and eigenfunctions of the Hamiltonians (j5J) and (|22j) 
particularly for highly excited Rydberg states we adopt the linear variational principle. Here 
the bound state solutions of the Schrodinger equation are expanded in a finite set of square 
integrable basis functions. Determining the expansion coefficients is equivalent to solving a 
generalized eigenvalue problem in case of non-orthogonal basis functions. The latter is done 
numerically by employing standard linear algebra techniques and routines. 

To accomplish the above we adopt spherical coordinates. The Hamiltonian (j5J) then 
becomes 

11 (3d 

H = — A r ibr sin 6 cos cj) I sin 2 9 cos9r— sin 3 9— 

2 r \ or 09 

b 2 b ( \ 

H — r 4 sin 4 9 sin 2 </> cos 2 cj) H — rsin# . (23) 

2 2 1 e -i4> o I 

With an additional Ioffe applied we have to consider the Hamiltonian ()22|) which reads in 
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spherical coordinates 



Hl = H-i^4j + —^ sin 2 o + El az . (24) 
2 86 8 2 v ; 



We utilize a Sturmian basis set of the form 

\nJ,m,s) = R^ k \r)Yr(9,6)\s). (25) 

These functions form a complete set in real and spin space but are not orthogonal. The 
angular part is covered by the well-known spherical harmonics Yf l {6,6) whereas the two 
spinor components are addressed by the spin-orbitals \s) — ||) or \[). For the radial part 
we employ 



RF\r) = ^^_^e-^(Cr)*L?(Cr). (26) 

with L^ k (r) being the associated Laguerre polynomials. The parameters k and ( can be 
adapted in order to gain an optimal convergence behavior in any spectral region. In par- 
ticular the non-linear variational parameter ( has to adapted such that it corresponds to 
the inverse of the characteristic length scale of the desired wavef unctions. Similar basis sets 
have been employed previously by several other authors ^| . 

The general expansion of an energy eigenstate \E) in a finite set of basis functions (|23j) 
reads 

\ E ) = y~] °nims \n,l,m,s) . (27) 

nlms 

From our knowledge of the symmetries of the system we can further specify the appearance 
of the expansion. In section IIHI we chose H , Y, z and P y P z I X yS2 to be the set of commut- 
ing operators whose eigenfunctions we want to construct. We now demand \E) to be an 
eigenstate of P y P z I xy S2- Exploiting the relations 

P y P z I xy S 2 Yr\]) = -ie-^ m (-l) l Yr\1) (28) 
P y PJ xy S 2 Yr\i) = e-^ m {-l) l Yr\[). (29) 
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we construct the following expansions for the four K-subspaces 

\E, +1) = ^Rni^nlmY^l+t 1 + ^nlmX^l^^ ID + R-n {CnlrnY^l+l 2 + d n i m Y 2 4m+4 ) IX)] (30) 

nlm 

\E, -1) = [Rn(a nlm Y 2 4 r +1 + b nlm Y 2 tt 3 ) IT) + Rn(cnimY 2 4 r +2 + d nlrn Y 4 ™+ 4 ) ID] (31) 

nlm 

\E, +i) = [Rn(a nlm Y 2 4 r +2 + b nlm Y 4 ™+ 4 ) IT) + RnicnimY^ 1 + d nlm Y 2 4 r +3 ) II)] (32) 

nlm 

\E, —i) = \Rn{o-nimY 2 4 l^i 2 + b n i rn Y 2l m+4 ) ||) + R n (c n i m Y 2 4m+1 + d n i m Y 2 4 ^ 3 ) ||)](33) 

nlm 

The eigenfunctions (I30H33P are automatically also eigenfunctions to T, z (see eq. (JHj) and 
()15))). Due to the structure of the spherical harmonics Y™ one has to ensure that \m\ < I. 
In our calculations the sums run over all valid combinations of n < N, I < L and m < M 
where the maximum indices N, L and M can be fixed individually. The expansion becomes 
exact if M, N, L — > oo. 

Performing the linear variational principle with one of the above expansions leads to a 
generalized eigenvalue problem Hv = ESv, where H and S are the corresponding matrix 
representation of the Hamiltonian ()23|) and the overlap matrix, respectively: 

H = (E,k\H\E,k) S=(E,k\ E,k). (34) 

The vector v contains the expansion coefficients a n [ m , b n [ m , c n \ m and d n i m . 

Due to the particular choice of the basis functions (J25|) the matrices H and S become 
extremely sparse occupied (S is a penta-banded matrix). In order to solve the generalized 
eigenvalue equation we utilize the so-called Arnoldi method together with the shift-and- 
invert method. We adopt routines from the ARPACK package. A more detailed description 
can be found in 

VII. RESULTS AND DISCUSSION 

In this section we analyze our computational results i.e. the eigenvalues and eigenfunc- 
tions obtained via the numerical approach described in the previous section. We discuss 
the spectra and expectation values of several observables as well as the properties of the 
electronic spin. Furthermore selection rules for electric dipole transitions as well as their 
strengths are derived. Results for the case of the additional presence of a homogeneous bias 
field are presented as well. 
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A. Spectral Properties 

With respect to the spectral behavior one can distinguish three regimes: the weak, the 
intermediate and the strong gradient regime each of which reveals individual characteris- 
tics. The appearance of these regimes is not determined by the gradient and the degree of 
excitation, i.e. energy, but by the scaled energy (see discussion in Sec. |TJ). For simplicity we 
will refer to the gradient as the relevant quantity characterizing the different regimes. All 
figures in this subsection show energy levels for manifolds belonging to rather small values 
for n (typically (n = 5 — 7)) and for large gradients (we cover the range b = 10~ 7 — 10~ 4 ) 
that are not accessible in the laboratory. This was done for reasons of illustration: Our 
observations and results equally hold for weaker gradients and higher n-manifolds (gradi- 
ents achievable for tight traps on atom chips are of the order of b = 10~ 8 ) which however, 
due to the high level density, are less suited for a graphical presentation. In the weak gra- 
dient regime the spectral behaviour is determined by the linear Zeeman terms. Although 
the principal quantum number n is not a good quantum number any given level can be as- 
signed to a certain n-multiplet. The levels split symmetrically around the zero- field-energy 
exhibiting the expected linear dependence on b. In figure |2k this is exemplarily shown for 
the n = 5-multiplet. 




FIG. 2: a: Splitting of the energy levels belonging to the n = 5 multiplet (k = ±l-subspace) with 
increasing gradient. The level structure is dominated by the linear Zeeman term. The splitting 
is linear and symmetric around the energy for b = 0. b: Intra n-manifold mixing of the n = 5 
multiplet in the k = 1-subspace. Due to the increasing dominance of the diamagnetic term the 
level splitting becomes non- linear. 
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The intermediate regime is characterized by the occurence of intra n-manifold mixing. 
Although neighboring n-manifolds are still distinguishable the levels now aquire a a nonlinear 
6-dependence which is due to the increasing importance of the diamagnetic term. Sub-levels 
belonging to different angular momenta mix and thus avoided level-crossings appear. The 
onset of this intermediate regime scales according to b oc n~ 6 . Figure |2fc> shows the regime 
of intermediate gradients of the n = 5-multiplet. Interestingly we observe here that this 
nonlinear behaviour in the /—mixing regime is very weakly pronounced for the atom in the 
side guide compared to an atom in a homogeneous magnetic field |2J. As we enter the 
strong gradient regime adjacent n-manifolds begin to overlap. The spectra are strongly i.e. 
nonperturbatively influenced by the diamagnetic term. Figure El shows this inter n-manifold 
mixing for the n = 6- and n = 7-multiplet where the strong coupling leads to large avoided 
crossings. The mixing threshold scales according to b oc n~^ (indicated by the dashed line 
in figure El). 



FIG. 3: Inter n-manifold mixing between the n = 6- and n = 7-multiplet in the k = 1-subspace. 
The mixing threshold is indicated by the dashed line. A large number of avoided crossings occur. 

B. Properties of the Electronic Spin 

1. S z expectation value 

In order to study the mutual influence of coordinate and spin space let us investigate 
the properties of the electronic spin. The x- and y-components of the spin operator obey 
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{S 2 , S x } = {S 2 , S y } = 0. Hence using (fTHj) we arrive at 

(S X ) = (Sy) = 0. 



(35) 



Only the expectation value of S z is non-zero in general. This is not obvious since the 
Hamiltonian (|2*3*j) does not contain an explicit dependence on S z . Figure 0^ shows the 
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FIG. 4: a: Expectation value of the z-component of the electronic spin operator for several excited 
states (b = 10 -7 ). b: Zoomed view of the n = 15-multiplet. The magnitude of (S z ) decreases for 
states possessing a large energy shift due to the external field. 

expectation value (S z ) for several excited states as a function of the principal quantum 
number n, which serves as an energetic label. The expectation values are arranged along 
vertical lines each of which belongs to a certain n-multiplet. With increasing degree of 
excitation these lines widen and begin to overlap as the inter n-mixing regime is reached. 
A zoomed view of the n = 15-multiplet is shown in figure EJd. We find states experiencing a 
large energy shift due to the external field thereby possessing a small S z expectation value. 
For the states shown in this figure (S z ) vanishes for n > 15.2 and n < 14.8. 
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2. Spatial Distributions of the Spin Polarization 



We now study the relative alignment of the electronic spin and the magnetic field. For a 
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(36) 



| M (r)| 2 + |d(r)| 2 

tFs_B(r) describes the spatial distribution of the spin polarization relative to the local mag- 
netic field. Wsb^ — 1 indicates the spin to be oriented parallel to the field whereas we find 
it antiparallel aligned for Wsb{t) — — 1- According to (JHrj|) Wsb^J) can be interpreted as the 
local expectation value of the cosine of the angle 7 between S and B. Since in a homogenous 
field the projection of the spin onto the field direction is conserved WsB{r) would be either 
+ 1 or —1 throughout the whole space. In the field of the side guide, however, we expect a 
much richer structure resulting from the coupling of the coordinate and the spin degrees of 
freedom. Figure El shows three tomographic cuts of a the spin polarization Wsb of the 83rd 
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FIG. 5: Tomographic cuts through the spin polarization Wsb of the 83rd excited state. The state 
belongs to the n = 8 multiplet within the k = 1-subspace (b = 10~ 7 ). The cuts are made at 
z = ±20 and z = 0. Positive and negative values are indicated white and black, respectively. We 
observe a rich pattern of different spin polarizations around the origin. From r ~ 60 on the nodal 
structure is replaced by a regular striped pattern varying periodically with the azimuthal angle </>. 

excited state in the k = 1-subspace. In the vicinity of the coordinate center we observe a 
large number of nodes. From p pa 60 on the complex nodal structure is replaced by a smooth 
regular pattern exhibiting a periodicity with respect to the azimuthal angle 4>. Here Wsb{^) 
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becomes almost independent of the z-coordinate. This feature seems to be induced mainly 
by the magnetic interaction which is invariant under translations along z. One identifies 
four sectors reminiscent of the quadrupolar structure of the magnetic field of the guides. In 
the present case we apparently have a anti-parallel alignment in the x — 0- and y = 0-plane 
and a parallel one between these planes. The densities are invariant under the operations 
P x P y and P z I X y which are equivalent to S z and P y P z I xy S2 when acting on real and scalar 
quantities. 

C. Electric Dipole Transitions 

We now consider electromagnetic transitions between electronic states in the framework 
of the dipole approximation. The transition amplitude between the initial state \i) and the 
final state |/) is then given by the squared modulus of the matrix element (i\D\f). In the 
length gauge D takes the forms D a ± = -^(x ±iy) = -^r sin 9e ±l< ^ and D w = z = r cos 9 for 
a - and 7r-transitions, respectively. 

Exploiting the symmetry properties of the P^P^J^^-eigenstates yields 

(E, k\ (PyP z I xy S 2 ) + z P y P z I xy S 2 \E', k') = kV (E, k\ z \E', k') = - (E, k\ z \E' , k') (37) 
which leads to the expression 

(W + 1) (E,k\z\E',k') = 0. (38) 

Here we have used (E, k\ {P y P z I xy S2) + = {k\k*. Apparently the matrix element for 7r- 
transitions can only be non-zero for the following combinations of K and k'\ 

tt : («, k') = (1, -1), (-1, 1), (z, -i), (-i, i) (39) 

The above shows also that the expectation value of the z— coordinate vanishes for any 
eigenstate i.e. we have (E, k\z\E, k) = For o" ± -transition one obtains in a similar way 

a + : («, k') = (i, 1), (1, -i), (-1, i), (-z, -1) (40) 
a' : («, «') = (-», 1), (1, i), (-1, -i), (i, -1). (41) 

Figure |H1 presents an overview of the allowed dipole transitions between the K-subspaces. 
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FIG. 6: Graphical representation of allowed dipole transitions between K-subspaces. The arrows 
point from k to k' . 
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FIG. 7: a: Dipole strenghts for 7r-transition from the ground state of the k = 1-subspace to 
excited states belonging to the k = —1-subspace (b = 10~ 7 ). The line for smallest A belongs to 
the n = 1 — > 7-transition. b: Zoomed view of the line belonging to the transition to the n = 12- 
multiplet. The line center is dominated by two sub-lines. The two bunches accompanying the line 
center at its left and right hand side possess a much smaller dipole strength. 

We have calculated the dipole strengths for transitions from the ground state to excited 
states. Figure U\ shows the results we obtain for 7r-transitions between the k — 1- and 
k! = —1-subspace. In figure UK we observe a general decrease of the dipole strengths with 
decreasing transition wavelengths. However, the decrease is not monotonous as it would be 
in the case of a homogeneous or a 3D-quadrupole field [r£|. One rather finds a modulation 
on top of the transition amplitudes where the n = 8-, n = 10- and n = 12-multiplet exhibit 
smaller dipole strengths than both of their neighbors. Figure [7jp shows a zoomed view of 
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the n = 1 — > 12 transition line. Its structure is dominated by two sub-lines located in 
the line center. The central bunch is almost symmetrically accompanied by two bunches 
of sub-lines located for smaller and larger wavelength, respectively. For cr + -transitions the 
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FIG. 8: a: Dipole strenghts for cr + -transition from the ground state of the k = — i-subspace to 
excited states belonging to the k = — 1-subspace (b = 10 -7 ). The line for smallest A belongs to 
the n = 1 — > 7-transition. b: Zoomed view of the line belonging to the transition to the n = 12- 
multiplet. The line consists of three bunches each of which consisting of a number of sub-lines. 
The line is dominated by two sub-lines one of each located in the left and right hand side bunch. 



dipole strengths are systematically decreasing with decreasing wavelength (figure|Ht). In the 
zoomed view (figure IHb) we also notice the structure consisting of three bunches of sub-lines. 
Again there are two dominating lines which are now located in the two outer bunches rather 
than in the central one. 



D. Magnetic Guide with a Ioffe field 

As discussed in section |TT] an additionally applied homogeneous field leads to severe 
changes of the symmetry properties of the atomic system. Apart from the lifting of the 
degeneracies also a significant influence on the electronic spin and the transition amplitudes 
have to be expected. 

Apparently there has to be a critical radius p c at which both fields are equal in strength. 
For a given gradient b and homogeneous field strength Bj it is given by p c = Taking 
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into account the scaling (p) oc n 2 we expect states with 




(42) 



to be equally affected by both fields. Hence, the states having n ■€ fi c or n > n c should 
be dominated by the homogeneous field or the field of the side guide, respectively. Figure 



FIG. 9: Expectation values of the z-component of the electronic spin at a finite homogeneous field 
strength (Bj = 10 -5 ) and a gradient of b = 10~ 7 . At low degree of excitation the homogeneous 
field dominates the electronic states. At this regime S z becomes an approximate constant of 



principal quantum number n c become increasingly dominated by the quadrupole field. As a result 
the expecation values tend towards (S z ) = 0. 

El shows the expectation values of S z for a gradient b = 10~ 7 and a homogeneous field 
strength Bj = 10 -5 . This yields the critical principal quantum number n c = 10. Indeed 
one finds for n <C 10 the expected dominance of the homogeneous field. In this regime (S z ) 
is approximately allowed to possess one of the two values ±|. This is due to the fact that 
S z becomes an approximate constant of motion. For n > 10 we observe the expectation 
values to move towards zero which is expected from the results shown in figure EJ We have 
to remark that since the symmetry S 2 persists the expectation values of S x and S y vanish 
even for finite strength of the homogeneous field. 

Apart from the spin expectation value also the spin polarization exhibits significant 
changes if a Ioffe field is switched on. For a sufficient high field strength or low degree 
of excitation (n < n c ), respectively, the structure of the electronic states is dominated by 




motion admitting (S z ) to possess one of the two possible values States lying above the critical 
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the Ioffe field. Here the spin is expected to be aligned with the homogeneous field. Since 
Wsb describes the projection of the electronic spin onto the direction of the side guide field 
which is perpendicular to the Ioffe field one expects Wsb to be approximately zero in this 
regime. Figure [TU] illustrates the SB polarization Wsb (equation (|3b|)) for the state shown 
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FIG. 10: Tomographic cuts through the spin polarization Wsb (equation of the 83rd excited 
state at a finite Ioffe field strength Bj = 10 -5 . The state belongs to the n = 8 multiplet inside the 
k = 1-subspace (b = 10 -7 ). The cuts are made at z = ±20 and z = 0. Positive and negative values 
are indicated by white and black, respectively. One observes large gray regions with Wsb ~ 0. 

in figure but for a Ioffe field strength Bj = 1CF 5 . The state is located inside the n = 8 
multiplet which lies below the critical quantum number n c = 10. Thus the states structure 
is dominated by the Ioffe field. As expected from the discussion above we observe large gray 
regions indicating Wsb = 0. The geometry of the side guide field is barely recognized for the 
cut made at z = 0. Unlike in figure El there are only small regions exhibiting a well-defined 
spin orientation that is dominated by the side guide, i.e either Wsb = _ 1 or Wsb = 1- 

Figure ITTk shows the dipole strengths for 7r-transitions from the ground state in the 
k = 1-subspace to various states in the k = —1-subspace. Compared to the Bj = case 
the dipole strengths are increased by approximately 70%. The transition strengths increase 
with increasing transition wavelengths. Again there seems to occur some kind of modulation 
as already seen in figure [?K but being less pronounced here. In the present case the n = 12- 
transition exhibits a larger transition amplitude than its neighbors. In figure ITTb we show 
a zoomed view of the line belonging to the n = 12-transition. Due to the presence of the 
homogeneous field a number of additional lines appear some of which are marked by an 
arrow. In contrast to the Bj = case the n = 12 line is dominated by a single sub-line 
originating from a transition induced by the presence of the homogeneous field. 
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FIG. 11: a: Dipole strengths for 7r-transition from the ground state of the k = 1-subspace to 
excited states belonging to the k = —1-subspace (b = 1CP 7 and Bj = 1CP 5 ). The line for smallest 
A belongs to the n = 1 — ► 7-transition. b: Zoomed view of the line belonging to the transition to 
the n = 12-multiplet. Several additional lines appear at finite homogeneous field strength (some are 
marked by an arrow) . The line center is dominated by a single sub-line emerging from a transition 
which is induced by the external homogeneous field. 

VIII. CONCLUSION AND OUTLOOK 



We have studied electronically excited hydrogen atoms located in a magnetic guide. In- 
cluding pseudo-potentials [l^ for the atomic core could be straight forwardly extended to 
describe e.g. alkali atoms. The magnetic guide represents a microtrap used to confine ultra- 
cold atomic systems. The motion of the valence electron has been described by an effective 
one-body approach. Both the coupling of the spatial degrees of freedom (para- and diamag- 
netism) as well as the spin degrees of freedom to the external field have been taken into 
account. The linear variational principle has been used to solve the stationary Schrodinger 
equation: Employing a Sturmian basis set enabled us to converge a large number of eigen- 
functions. 

A careful inspection of the Hamiltonian yields an amazingly large number of symme- 
tries involving both the spin and spatial degrees of freedom: We have found 15 symmetry 
operations of both unitary and anti-unitary character. This allows for a classification of 
the electronic eigenstates with respect to a complete set of commuting constants of motion. 
The latter involve the Hermitian S z -operator which is a combined spin and parity operator 
and the unitary but non-Hermitian operator P y P z I xy S2 which involves parity and permu- 
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tation operators. Employing specific anticommuting operators of this symmetry group we 
could prove the two-fold degeneracy of each energy level. This feature is indeed shown to 
be generic for spin-|-systems exhibiting certain symmetry properties. We have discussed 
how the symmetries are affected if an additional homogeneous magnetic field is applied in 
order to obtain a Ioffe-Pritchard type trap. In this case only 7 symmetry operations remain 
including E z , P y P z I xy S 2 and P x P z l xy S* 2 . 

Spectra have been investigated up to energies corresponding to a principal quantum 
number of n ~ 15. In the low gradient regime degenerate 72,-manifolds split up symmetrically 
around the zero field energy. For the intra-n-mixing regime only a very weak restructuring 
takes place inside any n-multiplet i.e. we observe only a minor nonlinear behaviour of the 
energies on the gradient. For even higher gradients the inter- n-mixing takes place where 
states belonging to adjacent multiplets begin to mix and avoided crossings dominate the 
spectrum. Scaling relations for both, the inter- and the intra-n-mixing have been provided. 

Effects due to the coupling of the spin and spatial degrees of freedom have been studied 
in detail. An analysis of the spin-field orientation has been performed by utilizing the 
distribution of the spin polarization. For electronic states in the magnetic guide Wsb reveals 
a rich nodal and island structure which is absent for an atom in a uniform field. Moreover 
an analysis of the S z expectation value has been performed. It has been shown that states 
being energetically strongly affected by the presence of the magnetic guide possess a small 
expectation value of S z . 

We have derived selection rules for the quantum number k belonging to the P y P z I xy S2 
symmetry operator for linear as well as circular polarized dipole transitions. Wave lengths 
and dipole strengths from the ground to Rydberg states were analyzed. In particular for n 
transitions we have found a global modulation of the transition amplitudes. The impact of 
the presence of an additional homogeneous magnetic field (along the wire involved in the 
set-up of the side guide) on several relevant quantities has been studied. This includes the 
S^-expectation values and the electric dipole transition amplitudes. 

Let us now comment on the approach chosen in the present work. Neglecting the fine 
and hyperfine structure of the atom as well as omitting the influence of the core scattering 
events represent, at least for certain species and regimes (high excitations !), certainly a 
good approximation to the true physical system. Another approximation is the fact that 
we centered the nucleus at the minimum of the field configuration. This is suggested by 
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our assumption that we have ultracold atoms with an extremely small kinetic cm. energy 
in tight traps leading to a well- localized atomic cm. Nevertheless, it is expected that the 
cm. motion blurs the effects ocurring for an atom with a fixed nucleus. Beyond this, 
it is well-known that already in the presence of a homogeneous magnetic field the cm. 
and electronic motions of atoms do not separate i.e. they perform an intimately coupled 



motion 
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24j . Then the immediate question arises how this coupling might 



look like in our inhomogeneous field configuration and in particular what its impact on the 
overall electronic motion is. To investigate this is a challenging task which needs careful 
consideration and clearly goes beyond the scope of the present work. 
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